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This work studies the feasibility of optimal control of high-fidelity quantum gates in a model of interacting two-level particles. One 
—I, particle (the qubit) serves as the quantum information processor, whose evolution is controlled by a time-dependent external field. The 
other particles are not directly controlled and serve as an effective environment, coupling to which is the source of decoherence. The 
7— H ■ control objective is to generate target one-qubit gates in the presence of strong environmentally-induced decoherence and under physically 
■ motivated restrictions on the control field. It is found that interactions among the environmental particles have a negligible effect on the 
If} ' gate fidelity and require no additional adjustment of the control field. Another interesting result is that optimally controlled quantum gates 
are remarkably robust to random variations in qubit-environment and inter-environment coupling strengths. These findings demonstrate 
the utility of optimal control for management of quantum-information systems in a very precise and specific manner, especially when 
the dynamics complexity is exacerbated by inherently uncertain environmental coupling. 
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1 Introduction 



The methods of optimal control are very useful for effectively managing various quantum systems pQ2| and 
are particularly important in situations requiring precise quantum operations, as is the case for quantum 
computation (QC) [3]. One of the most difficult problems of QC is that unavoidable coupling of the quantum 
' information processor (QIP) to the environment results in a loss of coherence. In recent years, significant 
attention was devoted to various methods of dynamic suppression of environmentally-induced decoherence 
in open quantum systems, including applications of pre-designed external fields (UEJE3E1E] and optimal 
control techniques pfMfTllfl^fT^fMfTH] . In a separate line of research, several works pBfl71fMfT9l[2nl[2T] 
considered the generation of optimally controlled unitary quantum gates in ideal situations where coupling 
to the environment can be neglected during the gate operation. 

The optimal control of quantum gates in the presence of decoherence still remains to be fully explored. 
Two recent works (22JE3] discussed specific techniques, involving optimizations over sets of controls op- 
erating in pre-designed "weak-decoherence" subspaces. We recently proposed [24] a different approach in 
which the full power of optimal control theory is used to generate the target gate transformation with the 
highest possible fidelity while simultaneously suppressing decoherence induced by coupling to a multipar- 
ticle environment. This method does not rely on any special pre-design of system parameters to weaken 
decoherence (e.g., using tunable inter-qubit couplings as in [22J or auxiliary qubits as in [23]); the only 
control used in our approach is a time-dependent external field. A similar application of optimal control 
was also recently considered in [25J for another model of a decohering environment. Optimization tech- 
niques were also applied recently to quantum error correction (QEC) [261127] . In contrast to QEC, our 
approach does not require ancilla qubits and is not limited to the weak decoherence regime. The optimal 
control of quantum gates can potentially be used in conjunction with QEC to achieve fault tolerance with 
an improved threshold. 
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In the previous work [24], we showed that optimal control fields, found by employing a combination 
of genetic and gradient algorithms, are able to produce high-fidelity quantum gates in the presence of 
strong decoherence. Optimal solutions revealed interesting control mechanisms that utilize dynamic Stark 
shifts to weaken coupling to the environment and control-induced revivals to restore coherence. In the 
present work, we extend the analysis of optimally controlled quantum gates to situations where (i) the 
environmental particles interact with each other and (ii) couplings between the QIP and environment 
and within the environment itself have randomly varied strengths. Taking into account these additional 
environmental effects makes our model more closely related to realistic quantum information systems (in 
particular, spin-based solid-state realizations of QC [28, 29, 301ET1E2] ) . An interesting finding is that the 
effect of inter-environment interactions on the fidelity of optimally controlled quantum gates is negligible 
(potentially allowing one to neglect certain interaction terms in the Hamiltonian). We also demonstrate 
that the optimal control generates quantum gates which are inherently robust to random variations in 
qubit-environment and inter-environment coupling strengths. 



2 Model system 

We use a model of interacting two-level particles (e.g., spin-half particles or two-level atoms), which 
are divided into the QIP, composed of one qubit, and an n-particle environment. The qubit is directly 
coupled to a time-dependent external control field, while the environment is not directly controlled and is 
managed only through its interaction with the qubit. The evolution of the composite system of the qubit 
and environment is treated in an exact quantum-mechanical manner, without either approximating the 
dynamics by a master equation or using a perturbative analysis based on the weak coupling assumption. 
The Hamiltonian for the controlled system, H = Hq + He + -ffint, nas the form (K = 1) 

n 

H = y^UjSjz - nC(t)S 0x - y^7ijSj • Sj. (1) 

i=0 i<j 

Here, i = labels the qubit and i = l,...,n label the environmental particles, Sj = (Si x , Si y , Si Z ) is 
the spin operator for the ith particle (Si = ^<Xj, in terms of the Pauli matrices), Hq is the sum over 
the free Hamiltonians ooiSi z for all n + 1 particles (uJi is the transition angular frequency for the ith 
particle), He specifies the coupling between the qubit and the time-dependent control field C(t) (/i is 
the dipole moment), and H mt represents the Heisenberg exchange interaction between the particles (7™ 
is the coupling parameter for the ith and jth particles). This model is particularly relevant to spin-based 
solid-state realizations of quantum gates (see, e.g., [28, 29, 30, 31, 32]). in which unwanted interactions exist 
due to impurities in semiconductor structures or usage of relatively dense lattices of spin-like qubits (e.g., 
electron spins in an array of quantum dots or electrons on liquid helium). 

In this work, we optimize one-qubit gates coupled to n-particle environments (n = 1, 2, 4, 6). For n = 2, 
the system can be modeled as a two-dimensional triangular lattice with the qubit qo coupled to two 
environmental particles e\ and &2'- 
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For n = 4, the system can be modeled as a two-dimensional lattice with the qubit qo at the center, coupled 
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Similarly, for n = 6, the system can be modeled as a three-dimensional lattice with the qubit at the 
center, coupled to six environmental particles. Since the evolution of the composite system is numerically 
exact and our optimization procedure (described in Section [5] and [23]) is iterative, we limit the number 
of environmental particles to n < 6 for computationally tractable simulations. 
We will first consider the case of well-specified coupling parameters given by 

J 7, for i = and j = 1, . . . ,n, , 
^ 1 j', for i = 1, . . . ,n — 1, and j > i, 



which implies that the qubit interacts with each environmental particle with the coupling parameter 7 and 
environmental particles interact with each other with the coupling parameter 7' (compare to [24] where 
7' = 0). Then we will consider a more general situation in which every coupling parameter 7^ takes a 
random value from a normal distribution. 



3 Distance measure for evolution operators 

Let U(t) G U(2 n+1 ) be the unitary time-evolution operator of the composite system and G G U(2) be 

the unitary target transformation for the quantum gate. The evolution is governed by the Schrodinger 

equation, U(t) = —iH(t)U(t), with the initial condition U(0) = I. The gate fidelity depends on the 

distance between the actual evolution U = f/(tf) at the final time if and the target transformation G. In 

order to achieve a perfect gate, it suffices for the time-evolution operator at t = if to be in a tensor-product 

form f7 op t = G ® <£, where $ G U(2 n ) is an arbitrary unitary transformation acting on the environment. 

Therefore, the following objective functional is proposed as the measure of the distance between U and 

G [33]: J = A n min||t/ — G ® $|| subject to $ G U(2 n ) (where || • || is a matrix norm on the space il#2»+i (C) 
'I' 

of 2 n+1 x 2" +1 complex matrices and A n is a normalization factor). Using the Frobenius norm, defined as 
||-A||pt = [trace (^4^)] 1/ ' 2 , and X n = 2~( n+2 )/ 2 , the distance measure becomes [33] 



J 



1 - 2A2trace ( y/QiQ 
2 



1/2 



rr'vv' 1 



(5) 
(6) 



r,r'=l 



where Q G M 2 ™ (C) and Q uu >, G rr ', and U rr i vv > are elements of the matrix representations of Q, G, and 
U, respectively. Since < J < 1, it is convenient to define the gate fidelity as F = 1 — J. An important 
property of this distance measure is its independence of the initial state. In contrast to some other distance 
measures [3J, J is evaluated directly from the evolution operator U, with no need to specify the initial state 
of the system. This property of J reflects our objective of generating a specified target transformation for 
whatever initial state, pure or mixed, direct-product or entangled. 
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4 Measure of decoherence and system parameters 

A useful measure of decoherence is the von Neumann entropy: S v ^(t) = —trace {p q (t) In [/9 q (i)]}, where 
p q (t) is the reduced density matrix for the qubit, p q (t) = trace en v (p(t))- For a pure state, 5" v n = 0, 
while for a maximally mixed state of a k- level system, SVn = ln(A;). The initial state used for the entropy 
calculations is |\£(0)) = [— )o ® ®r=i !+)«' where Si z \±}i = ±||±)j. The distance measure J of Eq. © 
is independent of the initial state and consequently so are the optimal control fields found for the target 
gates and the corresponding fidelities. Therefore, the specific choice of the initial state for the entropy 
calculations places no limitations whatsoever on the generality of the optimal control results. 

For the optimal control simulations below, the system parameters are chosen to ensure complex dy- 
namics and strong decoherence: values of j/lo are up to 0.02, which is significant for QC applications, 
and frequencies Ui are close (but not equal) to enhance the interaction. We define the unit of time and 
introduce a natural system of units by choosing the qubit frequency loq = 1 for all simulations (implying 
that one period of free evolution is 2n). The frequencies of the environmental particles are: u\ ~ 0.99841, 
uj 2 ~ 1.00159, u 3 » 0.96007, w 4 ~ 1.04159, u 5 » 0.87597, lu 6 « 1.14159 (see [M] for details). Imposing 
upper limits on the gate duration (tf < 60) and coupling parameter (7 < 0.02) places the dynamics of 
the uncontrolled system in the regime where decoherence increases monotonically with time (before the 
entropy reaches its maximum value of Svn ~ In 2). This dynamical regime approximates some of the effects 
that the QIP would experience from a larger environment, in particular, preventing restoration of coher- 
ence to the qubit by uncontrolled revivals. Thus, any increase in coherence may be attributed exclusively 
to the action of the control field. 



5 Optimization procedure 

Combined genetic and gradient algorithms are employed to minimize the distance measure J of Eq. ([5]) 
(or, equivalently, to maximize the fidelity F) with respect to the control field C(t). The target quantum 
logical transformation is the one-qubit Hadamard gate (an element of a universal set of logical operations 
for QC [3]). 

When a genetic algorithm is used, the gate fidelity F is maximized with respect to a parameterized control 
field C(t) = f(t) ^2 e Aicos (u>gt + 0^). Here, f[t) is an envelope function incorporating the field's spectral 
width and An, u>e, and 9i are the amplitude, central frequency, and relative phase of the £th component 
of the field, respectively. A combination of these parameters ("genes") represents an "individual" (whose 
"fitness" is the gate fidelity), and a collection of individuals constitutes a "population" (we use population 
sizes of ~ 250). 

Removing the constraints on the control field imposed by the parameterized form above provides the po- 
tential for more effective control of the system. In this case the optimal control field is found by minimizing 
the following functional [T7][l"8]: 

K = J + Re J ! trace { [U(t) + iH(t)U(t))B(t)} dt + | J ' \C(t)\ 2 dt. (7) 

Upon minimization of K, the first integral term constrains U{t) to obey the Schrodinger equation (B(t) is 
an operator Lagrange multiplier) and the second integral term penalizes the field fluence £ = J ' \C(t)\ 2 dt 
with a weight a > 0. Applying the calculus of variations to K with respect to B(t) and U(t) yields the 
Schrodinger equation for U(t) and the time-reversed Schrodinger equation for B(t): B{t) = iB(t)H(t), 
with an appropriate final time condition. The optimal field is found iteratively, using a gradient algorithm 
(see [21] for optimization details). 

Despite the lack of direct coupling of the control field to the environment, it can be shown that the 
composite system described by Eq. ([TJ is completely controllable (up to a global phase), as defined in [53] . 
However, the restrictions on the gate duration and on the shape of the control field limit the achievable 
fidelity. 
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Figure 1. The gate fidelity F versus the qubit-environment coupling parameter 7, for the Hadamard gates optimally controlled in the 
presence of ra-particle environments (n = 1, 2, 4, 6). Values of 7 range from to 0.02 in increments of 0.001. 



6 Results: Optimally controlled one-qubit gate with multiparticle environments 

Fidelities of one-qubit Hadamard gates optimized in the presence of n-particle environments (n = 1, 2, 4, 6) 
are presented in figure Q] for 7' = and various values of the qubit-environment coupling parameter 7. 
For a one-particle environment, the control fields optimized for the actual values of 7 result in fidelities at 
least above 0.9994. In particular, we obtain F > 1 — 10 -6 for 7 = (a closed system) and F ~ 0.9995 for 

7 = 0.02 (the strongest coupling considered). However, it becomes more difficult to counteract decoherence 
as the number of the environmental particles increases; as seen in figure [TJ for larger values of n the gate 
fidelity decreases more rapidly as 7 increases. We also find that applying a control field optimized for a 
one-particle environment (n = 1) to n > 2 results in a significant fidelity loss (up to ten percent of the 
original value). This demonstrates the dependence of the optimal control on the size of the environment. 

Optimal control field parameters, gate fidelity, and final-time entropy for the Hadamard gate coupled to 
n-particle environments (n = 1, 2, 4, 6, 7 = 0.02, and 7' = 0) are reported in table [TJ The fields are intense, 
with maximum amplitudes ranging from approximately 2.0 (for n = 1) to 3.9 (for n = 2). The exact time 
structure of the optimal field is not intuitive and is delicately adjusted to the particular control application. 
For example, control fields optimized for 7 = 0.02 are not only more intense than those optimized for 7 = 0, 
they also have very different structures. We also find that high-fidelity optimal solutions for 7 = are 
obtained for control pulse durations if < 12.0. In comparison, for 7 = 0.02, high-fidelity optimal solutions 
are obtained at longer pulse durations (cf. results reported in table [H e.g., if = 15.4 for n = 2 and 
if = 25.0 for n = 1 and n = 4). When the qubit-environment interaction is on, the control field is required 
to generate the target gate transformation and at the same time counteract decoherence. As described 
below, the control accomplishes the latter goal by restoring coherence to the QIP and therefore longer 
pulse durations are needed in the presence of environmental coupling. As mentioned above, there are also 
significant differences in the control fields optimized for different numbers of environmental particles. 

We also would like to explore how interactions among the environmental particles affect optimally 
controlled gate operations. Interestingly, we find that for a given n-particle environment (n > 2 and 
7 = 0.02), some optimal control solutions obtained for 7' 7^ (e.g., 7' = |7 = 0.0175) are essentially 
identical to the solutions found for 7' = 0. This means that no additional adjustment of the control field 
is necessary to account for the effect of inter-environment couplings. Moreover, applying these optimal 
control fields to the systems with 7' 7^ yields approximately the same gate fidelity, as seen for 7' = 
(the actual decreases in the fidelity observed for 7' = ^7 and 7' = ^7 are of the order of 10 -5 for n = 2 
and 10 -4 for n = 4). 

Control mechanisms can be better understood by examining the decoherence dynamics of the qubit. 
Figure [2] shows the time behavior of the von Neumann entropy of the qubit for the optimally controlled 
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Tabic 1. Optimal control field parameters (the maximum field amplitude 
A raax , control duration tf, and field fluencc £), gate fidelity F, and final-time 
entropy S v j^(t{) for the Hadamard onc-qubit gates coupled to n-particle en- 
vironments (n — 1,2,4,6 and 7 — 0.02). The initial state for the entropy 
computation is |^(0)). 
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Figure 2. The von Neumann entropy S v s(t) versus time, for the optimally controlled Hadamard gates coupled to two- and 
four-particle environments (n = 2,4, 7 = 0.02, and 7' = 0.0175). The initial state for the entropy computation is |^(0)). 



evolution, with 7 = 0.02, 7' = 0.0175, and n = 2,4. The difference between the entropy values for 
7' 7^ and 7' = is extremely small, implying that the same control mechanism works in both cases. 
We observe that the optimal control dramatically enhances coherence of the qubit system in comparison 
to the uncontrolled dynamics. Decoherence is suppressed by the control at all times, but especially at the 
end of the transformation. For example, for the Hadamard gate with 7 = 0.02, S v ^(tf) < 10~ 7 in the 
presence of a one-particle environment, which means that at t = tf the qubit and environment are almost 
uncoupled. While the value of S V N(if) increases with the environment size (see tabled]), the optimal control 
is still able to achieve a significant degree of decoupling at the final time. Inspecting eigenvalues of the 
controlled Hamiltonian, we find that the intense control field creates significant dynamic Stark shifts of 
the energy levels. This effect is mainly responsible for reducing the qubit-environment interaction during 
the control pulse. However, achieving extremely low final-time entropies and correspondingly high gate 
fidelities requires the employment of an induced coherence revival. In the uncontrolled system, revivals 
occur at times much longer than tf, so that the strong coherence revival observed at t = tf is induced 
exclusively by the control field. However, as the complexity of the composite system increases, it becomes 
more difficult to induce an almost perfect revival; therefore, the gate fidelity and final-time coherence 
decrease as n increases. Similar results are found for a two-qubit controlled-NOT gate optimized in the 
presence of a one-particle environment, in which case we observe |24j a longer pulse duration (tf = 121.1) 
and smaller fidelity {F ~ 0.9798) than for one-qubit gates. 
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Tabic 2. Fidelity and entropy data for the onc-qubit Hadamard gate, obtained when the 
control field optimized for a specified set of the system parameters is applied to an ensemble of 
systems with normally distributed variations in the coupling parameters 75* . Columns of F and 
SvN^tf) contain fidelity and final-time entropy values, respectively, obtained for the specified 
coupling strengths: 7 — 0.02 and 7 — c-y (c — 0, 1/2,7/8). Columns of F and SVn contain 
mean values of the fidelity and final-time entropy, respectively, over the ensemble, while ap- 
and &s vN are the respective standard deviations. 
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7 Robustness of optimally controlled gates to coupling-strength variations 

In realistic quantum systems, the strength of coupling between particles not always can be accurately 
measured and is also subject to fluctuations caused by noise or variations due to imperfect manufacturing 
or preparation. Therefore, it is important to explore the effect of general symmetry-breaking variations 
in the coupling parameters 7^ on operations of optimally controlled quantum gates. Given the one-qubit 
Hadamard gate as the target transformation and a fixed number n of environmental particles (n = 2,4), 
we find the optimal control field for a specified set of coupling parameters jij of Eq. with 7 = 0.02 
and 7' = 07 ( c = 0, ±,§). Then we apply this control field to an ensemble of systems with normally 
distributed variations in the coupling parameters 7^ and analyze how the uncertainty in the coupling 
strengths affects the gate fidelity and final-time entropy. Although the dependence of F and S v ^s(t{) on 
the coupling parameters is non-linear (which implies that the corresponding distributions of F and S^N^f) 
will not be normal), our statistical analysis employs only mean values and standard deviations, given by 
F = L~ l Ylp=i Fp an d <?F = [L^ 1 J2p=i(Fp ~ F) 2 ] 1 / 2 , respectively, for the gate fidelity F, and similarly 
for the final-time entropy S v -^(t{). The summation is over all elements of the ensemble (ensemble sizes L 
of the order of 10 5 are used in the calculations). 

For each element of the statistical ensemble, the value of each qubit-environment coupling parameter 
7oj (j = 1, . . . , n) is randomly selected from the normal distribution with a mean 7 = 0.02 and a standard 
deviation 07 = 7/8 = 0.0025. Analogously, for non-zero inter-environment coupling^] the value of each 
coupling parameter 7y (i = 1, . . . , n — 1, j > i) is also randomly selected from the normal distribution with 
a mean 7' = 07 and a standard deviation ay = j'/8 = C7/8. The statistical analysis of the corresponding 
fidelity and final-time entropy distributions is reported in table [2] (for n = 2,4), and frequency histograms 
of these distributions are shown in figure [3] (for n = 4). These results demonstrate a high degree of 
robustness of the optimally controlled gates to relatively large variations in the qubit-environment and 
inter-environment coupling strengths. For given values of n, 7, and 7', on average there is just a minuscule 
decrease in the fidelity and final-time coherence due to the coupling strength variations, and the relative 
width of the fidelity distribution, opjF , is smaller than 07/7 by several orders of magnitude. Based on 
the data in tabled we also observe that the standard deviations of and as vN increase with the number of 
the environmental particles. Table [2] also helps us to see that inter-environment couplings have very little 
effect on the gate performance, although their influence slightly increases with 7'. 



In the case of zero inter-environment coupling, c = 0, all zero values of 7ij are left unchanged. 
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Figure 3. Frequency histograms for (a) the gate fidelity distribution and (b) the final-time entropy distribution. These distributions 
are obtained when the control field optimized for the Hadamard gate with n = 4, 7 = 0.02, and 7' = (7/8)7 = 0.0175 is applied to an 
ensemble of systems with normally distributed variations in the coupling parameters ■jij, as explained in the text. Table [2] reports 
statistical data for these distributions. Note that the sub-plots have different scales of the axes. 



8 Conclusions 



This work demonstrates the importance of the optimal control theory for designing quantum gates, es- 
pecially in the presence of environmentally-induced decoherence. The model studied here represents a 
realistic system of interacting qubits with uncertain coupling strengths and is relevant for various physical 
implementations of QC. Very precise optimal solutions obtained in the presence of unwanted couplings 
reveal control mechanisms which employ fast and intense time-dependent fields to effectively suppress 
decoherence via dynamic Stark shifting and restore coherence via an induced revival. In addition, these 
optimal solutions exhibit a significant degree of inherent robustness to random variations in the coupling 
strengths. It is also found that optimally controlled gate operations are practically unaffected by interac- 
tions between the environmental particles. These results further support the use of laboratory closed-loop 
optimal controls in QC applications. 
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